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Abstract 

We introduce a Hamiltonian description for the dynamics of several rigid bodies 
in teracting with point vortice s in an inviscid, incompressible fluid. We adopt the idea 
of lVankerschaver et al.l |200^ to derive the Hamiltonian formulation via symplectic re- 
duction of a canonical Hamiltonian system on a principle fibre bundle. On the reduced 
phase space we determine the magnetic symplectic form directly, without resorting to 
the machinery of mechanical connections on principle fibre bundles. We derive the 
equations of motion for the general case, and also for the special Lie-Poisson case 
of a single rigid body and zero total vorticity. Finally we give a partly degenerate 
Lagrangian formulation for the system. 



1 Introduction 

The Hamiltonian dynamics of a single rigid body of arbitrary shape in a 2-dimensional 
in viscid and incompre ssible fluid interacting with n poi nt vortices has first b een formulated 
by IShashikanth ( 20051 ). The system was also studi ed by Borisov et al. ( 2007 ). dropping the 
restriction of zero circulation around the cylinder. 'Vanker schaver et al. ( 2009 ^) showed that 
both formulations are indeed equivalent: Both can be derived when regarding the dynamics 
as geodesies on a Riemannian manifold, in the spirit of Arnold's geometric description 
of fluid dynamics. The manifold here is the cartesian product of SE(2) with a subset 
of volume-preserving embeddings of the initial fluid conflguration into (compatible with 
the time-dependent pose of the body), with kinetic energy as the Riemannian metric. When 
reducing the system to fluid velocity flelds generated by point vortices, one obtains a flnite- 
dimensional phase space with magnetic symplectic form, which yields the coupling between 
rigid body and point vortex motion. 

In the absence of vortic i ty, th e dynamics of several rigid bodies in a fluid has been 
studied b y Nair and Kansol (|2007l ). Their work is an extension of the original work by 
Kirchhoflj(|l87fl[ ). who discovered that the kinetic energy of a surrounding potential flow can 
be incorporated into the kinetic energy of the rigid motion as added mass. In contrast to the 
case of a single rigid body, the kinetic energy of the potential flow around several moving 
obstacles is no longer a constant quadratic form, but depen ds on the relativ e pose s of the 
different bodies. Using this kinetic energy as a Lagrangian, iNair and Kansol (j2007r ) derive 
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the equations of motion for a number of configurations, and study them analytically and 
numerically. 

In this paper we introduce the Hamiltonian dynamics of several rigid bodies interacting 
with point vortices, for the case of zero circulation around the individual bodies, but arbi- 
trary strengths of the point vortices. This system has been known only f or the case of a single 
rigid b ody. We adopt the description of the reduced phase space from IVankerschaver et al. 
( 2009t ) and extend it to the case of several rigid bodies. On the reduced phase space we de- 
termine the magnetic symplectic form directly, using only general properties of the magnetic 
symplectic form, and Kelvin's circulation theorem. From there we derive the equations of 
motion for the general case, and also for the special Lie-Poisson case of a single rigid body 
and zero total vorticity. Finally we give a Lagrangian description of the system. 

2 Physical Model 

2.1 Rigid Body Configuration 

Assume that we have a collection of n rigid bodies B = {Bi, . . . , Bn\ moving in an inviscid, 
incompressible fluid occupying the whole of . The bodies move rigidly, so their configura- 
tion is described by a curve 1 1— )> 5 in SE(2)", and B is the image of the initial configuration 
B at i = 0: B = g{B) = {51(^1), . . . ,g„(B„)}. We denote the body velocity by rj e se(2)", 
i.e., g = gfl. For the angular and linear velocity components of flj we write (ujj^Vj), i.e., 
— ujj X z + Vj. 

2.2 Fluid Configuration 

The time-dependent fluid domain J- (including body boundaries) is covered with a fluid 
which is at rest at infinity and whose motion is given by the time-dependent fluid velocity 
field u. Its vorticity field lu = curlu is zero everywhere, except for isolated point vortices 
7 = {71, . . . ,7m}. There the vorticity field is concentrated in a delta-function-like manner. 
The circulation around each vortex is constant in time (due to Kelvin's circulation theorem) 
and measures the strength Ti of the vortex 7^. We assume zero circulation around the 
individual bodies, and impose no-through boundary conditions. I.e., the normal component 
of the velocity field must coincide with the body boundary normal velocity, while the tangent 
velocity is arbitrary: 

{u{z),nj{z)) = {gjgj\z),nj{z)), for z G dBj. (1) 
2.2.1 Hodge-Helmholtz Decomposition 

We will now construct the fluid velocity field u for a given configuration {g,g,j). In the 
absence of bodies, the fluid velocity field whose vorticity is given by the point vortices 7 
with strengths F is determined by the Biot-Savart law: 

This defines a vector field on the whole of E^, and we denote its restriction to the fluid 
domain T hy u-y = ubs\t- 
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When bodies are present, the velocity field u-, makes fluid particles move across the 
body boundaries. To fix this, we first construct a potential field ui — grade/)/ on !F which 
compensates the normal flux of on the body boundaries, thus satisfying the boundary 
condition 

(u/(z), nj(z)) = — (u-y(z), nj(z)), for z £ dBj. (3) 

The subscript / reflects the fact t hat ui can be re presented as image vorticity inside of the 
bodies or on their boundaries, see ISafFmanI (jl992l §2.4). The potential 0/ of uj is uniquely 
determined by the Neumann problem 

^(z) = -{uJz),n(z)), for z G dT, A(t)i(z) = 0, lim grad0/(z) = 0. (4) 
an z->oo 

Now U|| — Uj + uj satisfies the boundary condition (z), n(z)) = on dJ-. In other words, 
it is the correct fluid velocity field as long as the bodies are at rest. When they move we 
achieve boundary condition ^ by adding the potential field = gradc/ig, obtained from 
the Neumann problem 



OB 



dn 



{gj{z)) ^ {nj{z),nj{z)), ioT z edBj, A(/)b(z) = 0, lim grad0B(z) = 0. (5) 



The superposition u = u^ + ui + us is the unique fluid velocity fleld which satisfies boundary 
condition ([T|), has zero circulation around the individual bodies, vanishes at infinity, and 
whose vorticity field is given by the point vortices 7^ with strengths F^. 



2.2.2 Stream functions 

Equivalently we can represent the velocity fields u/ and in terms of their stream functions 
ipi and ipjs- I-G-j u = grad0 = Jgradi/i. When we determine the symplectic form later on, 
we require the falloff behaviour of tpi, tpB and the variations dg^pj, dgipis with respect to the 
body configuration g. 

Lemma 1. The functions ipi, ipB, 9gipi and dgif^B ^{\^\^^)- 

Proof. Let u be uj or ug, and ip the stream function of u. We represent -0 as a single layer 
potential with density r, i.e., 

V'^(z) = * T{y)log\\y - gT^{z)\\dy, V(^) = V . V'^ 

JdB, ■' ^3 

Note that we can also view r as the strength of a vortex sheet (i.e., a distribution of point 
vortices) on dB, which generates u via the Biot-Savart law ([2]). Since the circulation of 
u = gradt/) around any Bj is zero, the total density of t on each dBj, as well as its variation 
with respect to g, vanish: 

f T{y)dy=(b {dgT){y)dy ^0. (6) 

J dBj JdBj 

This implies ip-'{z) — ipiz) = C'(|z|^"'^). For the variation with respect to g we obtain 
dgr{z) = ((f Tiy) J ~ dy,-dggT\z)) + I {dgT){y)\og\\y ~ g-\z)\\ dy. 
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The second term is also a single layer potential and Odz]"^) because of ©■ The boundary 
integral in the first term is the velocity field generated by the image vorticity on dBj , rotated 
by — 7r/2. Thus it is 0{\z\~'^) because of (O, and so the first term is as well. □ 

2.3 Kinetic Energy 

The kinetic energy of the system is the sum of kinetic energy associated to the rigid body 
motion and the kinetic energy of the fluid velocity field. We will stick with the common 
practice to only consider the finite part of the kinetic energy of the fluid, ignoring the infinite 
self-energy of point vortices. 

The kinetic energy of the motion of a single rigid body is given by its mass-inertia 
tensor A^^: Tm, — ^j)- Here we have identified ilj with a 3-vector containing 

Sj's angular and linear velocity. Summation over all bodies gives the kinetic energy of the 
motion in SE(2)" as a quadratic form on the Lie algebra se(2)": = ^{T,j\4T), where 
Ai contains the individual mass-inertia tensors on the diagonal. By left translation, this 
defines a left-invariant metric on SE(2)". 

The kinetic energy of the fluid velocity field u is 

U \wr = \j I \\ub\\\ (7) 

since u\\ and L^-orthogonal. The kinetic energy associated to can be expressed in 
terms of added mass as a quadratic form on the body configuration space SE(2)". We 
denote its representing matrix, the added mass tensor, by ICg: 

1 r „ 1 



j^\\usr = \{n,iCgn). 



For a single rigid body this g oes back to Khchhofj (l l87Cll ). while the case of more than 



one body was worked out by iNair and Kansol ( 2007 ). Together with the kinetic energy 



associated to the rigid body motion, we define the Kirchhoff tensor K — M + ICg. The 
corresponding energy, 

can be viewed as a function on the cotangent bundle r*SE(2)": Define P — Kil, then 

H,s = l{K-'P,P) (8) 

is the Hamiltonian for the dynamics of rigid bodies in potential flow. 

The kinetic energy associated to M|| contains infinite self-energy terms that we ignore. 



The finite part is given as the negative of the Kirchhoff- Routh function Wg ()Linl . 11941 
Shashikanthl . l2005l) : 



WG = lY.^^UI{l^)+J2i'M] ^ (9) 

where ipj^^z) — —Tj log \\^j — z|| is the stream function of the point vortex 7-,-. Note that 

^^^^G = -^^JU\^^(^^)■ (10) 

The total kinetic energy is of the coupled system is 

B^Bb-Wg^Tb- Wg- (11) 
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3 Hamiltonian Formulation 



In analogy to Arnold's geometric description of fluid dynamics (jArnoldl . 119661 ) . the dynam- 
ics of rigid bodies interacting with a surrounding incompressible fluid can be viewed as 
a geodesic problem on a Riemannian manifold. The kinetic energy defines a Riemannian 
metric on the configuration space, and the geodesies satisfy Hamilton's equations for the 
kinetic energy as the Hamiltonian, and the cota ngent bundle with canonica l symplectic 



structure as phase space. This insight is due to IVankerschaver et al. (I2OO9I) . From the 



general cot angent bundle descrip tion, the authors use the framework of cotangent bundle 



reduction ([Marsden et al.l . 120071 ) to obtain a reduced Hamiltonian system with magnetic 



symplectic form for the case of a single body in a fluid whose vorticity field is concentrated 
at point vortices. 

We will make use of the fact that the reduction by IVankerschaver et al. I (I2OO9I) does not 
rely on the fact that only a single rigid body was considered. However, when the authors 
derive the precise form of the magnetic symplectic form, they make use of the assumption of a 



single (and even circular) rigid body. At this point we will deviate from IVankerschaver et al 



and determine the magnetic symplectic form directly. In contrast, they compute the 
magnetic symplectic form using the machinery of mechanical connections (the Neumann 
connection here) on principle fibre bundles. 

The Hamiltonian formulation of Vankerschaver et al. I (I2OO9I) is an extension of Arnold's 
original work (jArnoldl . 119661 ) , where he described the motion of an incompressible inviscid 
fluid in some fixed fiuid domain as a geodesic on the group Diff vol (J^) of volume-preserving 
diffeomorphisms on J-. In the case when the fluid interacts wit h rigid bodies, the fluid 
domain is no longer flxed. The idea of Vankerschaver et al. ( 20091) is to consider the space 
Embvoi (-7^'' , ) of volume-preserving embeddings of an initial reference conflguration J-^ 
into instead of Diff vol (-7^) • Any incompressible fluid motion is then described by a curve 
/: t (pt which satisfies the condition ipf{J''^) = 7^*, where 7^* denotes the time dependent 
fiuid domain at time t. The configuration space of the coupled system isQ = SE(2)" x (Q-^ C 
EmbvoK-T^*^, ]R^)), and the dynamics is a canonical Hamiltonian system with kinetic energy 
as the Hamiltonian and T*Q with canonical symplectic structure as phase space. 

The kinetic energy is invariant under volume-preserving diffeomorphisms of the initial 
fluid configuration J^'^ {particle relabelling symmetry), i.e., the symmetry group Diff vol (7^°) 
acts from the right on Q'^, and thus on Q. This action makes Q into the total space of a 
principle fibre bundle over SE(2)". This principle fibre bundle structure allows to follow the 
famous Kaluza-Klein approach to determine the Hamiltonian dynamics. In order to factor 
out the Diffvoi(7^'^)-symmetry one needs to fix a value of the associated momentum map. 
This is where the assumption is used that vorticity is concentrated at m point vo rtices. The 



Vankerschaver et al. I (l2009l §4.2), and 



reduced phase space is Qred = T*SE{2Y x E^™, see 
the dynamics is given using the reduced Hamiltonian together with a reduced symplectic 
form a red on Qred- The reduced symplectic form Ured is the sum of the canonical symplectic 



Vankerschaver et al. I (I2OO9I 



form on T*SE(2)" and a magnetic term ap on SE(2)" x M"^", see 
Eqn. (44)). 

Be fore we determine we summarise the results that we adopt from lVankerschaver et al 
(|2009[ ): 



Theorem 1. The dynamics of n rigid bodies interacting with m isolated point vortices is a 
Hamiltonian system. The phase space of the system is 



Q 



red 



T*SE(2)' 
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where T*SE(2)" is the cotangent bundle of the rigid body configuration, andM.^"^ is the phase 
space for m point vortices. A point {p,g,^) € Qred consists of the rigid body momentum 
p S T*SE(2)", the rigid body configuration g G SE(2)", and the point vortex state 7. The 
momentum p = Ad*-i P encodes the rigid body velocity g = Adgfl (using the Adjoint 
representation o/SE(2)"J through P = K~^Q,. The Hamiltonian of the system is the kinetic 
energy as a function on phase space: 

H = -Wg + Hb. 
The symplectic form Gred for the system is 

^red = CTcan + TTpCT^, 

where Gcan is the canonical symplectic form on the cotangent bundle r*SE(2)", and TT*ap 
is a magnetic term; I.e., it is a two-form on SE{2)" x M^™, lifted to Qred through the 
canonical projection tTj, : (p, g, 7) 1— > {g, 7). 

Proof. This has been proven bv IVankerschaver et all ( 20091 §4). We emphasise here that 



the proofs do not rely of the fact that only a single rigid body was considered. □ 
It is now where deviate from lVankerschaver et aL 1 (I2OO9I) in order to obtain the magnetic 



term for the case of several rigid bodies. While they compute directly from the curvature 
of the Neumann connection, we determine it form from general properties of a magnetic 
symplectic form, and the limiting behaviour of the system for large distance between bodies 
and point vortices. Together with Kelvin's circulation theorem, this allows to determine the 
magnetic symplectic form directly from first principles, without resorting to the machinery 
of mechanical connections on principle fibre bundles. 

We will now study the behaviour of the coupled dynamical system in the limit when 
rigid bodies and point vortices are far apart. Assume that both bodies and point vortices 
are contained in two disjoint disks of finite radius, and let d denote the minimal distance 
between these disks. We will now show that in the limit d ^ 00 the system decouples: The 
rigid bodies move as if no point vort ices were present (so-called hydrodynamically coupled 



rigid bodies) (jNair and KansoL 120071 ) . and the point vortices move as if no boundaries were 
present. Hydrodynamically coupled rigid bodies are a canonical Hamiltonian system on the 
cotangent bundle T*SE(2)", with Hamiltonian Hb, see Equation (jS]). The dynamics of m 
point vortices 7^ (with strengths Ti) in the plane is a Hamiltonian system with phase space 
M^™, and symplectic form 

c^7(7'7) = y^.r* det(7,7). 

The phase space R^™ is a coadjoint orbit of the initial vorticity field under the action of 
the formal Lie-group Diff vol (K^ ) of volume-preserving diffeomorphisms of R^, and is the 
Kirillov-Kostant-Souriau form on that orbit. The Hamiltonian here is 

The decoupling can be expressed as follows: 

Lemma 2. The dynamical system decouples in the limit d — S- 00, i.e., 

lim ap = n*a^, (12) 

where TTg is the canonical projection (3,7) 7. 
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Proof. We lift the two uncoupled systems to Qred through canonical projections, and sup- 
press this lift in our notation. Now consider the difference of Hamilton's equations of the 
coupled and the two uncoupled systems. We have 

(Jred - (Jean - = ap - (J^ and H - Hb - = -Wq ~ ^7 ~^ X! ^ii/'/lTi)) 

i 

and obtain 

c^/3((5,7), (3,7)) - ^^7(7: 7) = d(i/ -Hb- Hy)Cg,j) 



The right hand side vanishes the limit d — >■ 00 because of Lemma [T] and the construction of 
ui, see the Equation (jH). □ 

Using the fact that point vortices are frozen into the fluid flow due to Kelvin's circulation 
theorem, i.e., ji — u{-^i), we can now determine a primitive one-form for ap: 

Theorem 2. The one- form 

i i 

is a primitive one-form of the magnetic term ap of the reduced symplectic form Ured on 
Qred- The components of^i^ are basis stream functions ofu^, i.e., tpeiz) ~ 

Proof. We can write cTred as 

CTred((p,ff,7),(P,5,7)) = f^can((p,5),(p,g)) + ( , ^ (^^^ ) ^ (13) 

" v ' 

where Mi, and depend only on g and 7, and Jb and J-y are screw-symmetric matrices. 
Now consider a solution {p{t), g{t),^{t)) to Hamilton's equations. In particular we have 
fT.ed((p,5,7), (0,0,7)) =di/(0, 0,7), i.e., 

(7,-A/*-g-f J^-7) =^^r,(J^.||(7,),7.)- (14) 

Neither ap nor the right hand side depend on p, so the above equation holds for any solution 
in the p-fibre over (17, 7). So we can consider p = which implies 5 = and thus ji — U|| (7;). 
This determines J^: the i-th diagonal block is TiJ. Plugging into ((Ti)) gives 

(7,-M| •<?) = - ^r,(Jue(7,),7.) =Er^(g^-ad(*B(^»)'")'^») 

i i 

= {n,{d,p^)-^)^{g,d,{Ad;^.p^)-j), ^'^^ 

-Ml 
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which determines Mj. From J-^ and M/ we can find a primitive one-form of as follows. 
Look at the restriction of ared to the R^™-factor. Since does not depend on g we have 

a,ed((0, 0, 7), (0, 0, 7)) = <p,3)'^7(7> 7), 

with (T-y the Kirillov-Kostant-Souriau form on R^™ and 7r(p : (p, g, 7) 1— 7 the canonical 
projection. Since cr^ is an exact two- form on M^™, we may assume that up possesses a 
primitive one-form of the form 0/3(7,5) = {-^{911)^9) + (-^(7)17) ^-nd compute 

da((<?, 7), (ff, 7)) = (g, [(a^AT^-ag^ g - d^A 7) (7, P^B^^a^ 7 + g). 
Equation ([T5|) determines A up to a contribution A that is independent of 7: A = P'~^ + A, 



c _ A^* _ pC ^^^^^^ 



with P'^ = Ad*-i P'^. Then 



= (a^P^)* - a^P^ + {dgAf - dgA, 

and, when considering the limit d ^ 00 using the Lemmas [1] and [2l we get 

lim ({dgAf - dgA) = 0. 

It follows that A does not contribute to at all because it is independent of d, and we may 
choose A = 0. It remains to check that 07(7) is a primitive one-form of cTj on M^"*, which 
completes the proof. □ 

3.1 Equations of Motion 

We can now derive the equations of motion by extremizing the action 
rti rti 



f ' ared -Hdt^ f\p,n)+a^ (7) + (P^,n) +Wg - He dt, 

J to J to '- v ' 

Wb 



(16) 



where Ored is a primitive one- form of ared- Note that W = Wb + Wq is t he Ki rchhofF-Routh 
function extended to rigidly moving boundaries, see IShashikanth et al. I (I2OO2I) . 

Let P be an arbitrary variation of P, while g and 7 are variations with fixed end points. 
We obtain 

d^(P, 0, 7) = f\p,n- K-'P) + (77(7, 7) + d^W{^) dt, 

Jto 

which just restates that P = Kfl and •ji = u{ji). Now consider g = gA, which implies the 

o 

Lin-constraint Q — ad^ A + A: 
rti 



d5(0,5,0)= / \n,P + P^)+dgWC9)-dgHBC9)dt 

Jto 

= r (A, -(P + pC-) + ad*,(P + pC) + Ad; dg{W - Hb)) dt. 

J to 



We obtain the following 



8 



Corollary 1. The equations of motion for n rigid bodies interacting with m point vortices 
in an inviscid, incompressible fluid are: 

Ph = Kn + p^, 

Ph = adi,PH + Ad; dg{w~HB), ^^^^ 

g = gn. 

3.1.1 The Lie-Poisson case 

A particularly interesting special case is the one of a single rigid body with point vortex 
strengths that sum up to zero. In this case the equat i ons of motion are a L ie-Poisson system. 
This has been proven earlier in ( Shashikanthl . 120051 : Borisov et al. , 2007 ) , and we will now 



derive it from our general equations of motion (I17p . by expressing them in the body frame of 
reference. This amounts to factoring out the SE(2)-symmetry of the system. Note that this 
symmetry also exists in the case of several bodies, and one could factor out the symmetry 
by expressing the equations of motion in the frame of reference of one particular body. 

In the case of a single body we can express the Kirchhoff-Routh function in the body 
frame of reference: W{li, . . . , Im) = W^(7i, ■ • ■ , 7m), where li = g~^(7i) are the point vortex 
positions in the body frame of reference. Here we have 

duW = -r,Ju{k), 

where u is the fluid velocity field u expressed in the body frame of reference. Now it is 
important to notice that all the stream functions involved in W are independent of g. Thus, 
together with Lemma [31 we obtain the following for the variation of W with respect to g: 

d,W{gA) = r,{Ju{h), Mh)) - E ( ^ ^ ) ' a) = i-Pi + ad?, Pi, A), 



with 



\hr/2 
Jk 



The Kirchhoff tensor does not depend on g in the case of a single rigid body, so the equations 
of motion are 



PH^&di,PH, k^^n{h) + n{k), Ph = Kn + p'^ + Pi. (is) 



The sum P'-' + Pi is the momentum of the fluid velocity field, see Ishashikanth ( 20051 ). and 
the equations of motion agree with the ones derived therein. 

Lemma 3. Let k and u denote point vortices and the fluid velocity field, both in the body 
frame of reference. Under the assumption that the vortex strengths sum up to zero, the 
following equality holds: 

Proof. Compute Pi, using li — —^{k) + u{li) and J2i^i — 0- ^ 
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4 Lagrangian Formulation 

We obtain a Lagrangian formulation for the dynamical system by considering the action (|16p 
as a fmictional on the tangent bundle of the configuration space SE(2)" x R^™: 

ti 

a^{'^) + W + Ti3dt. (19) 

When extremizing S we obtain the same equations of motion, since dgHg = —dgT^. In this 
way we have proven the following 

Corollary 2. The coupled system of rigid bodies and point vortices is a Lagrangian system 
with configuration space SE(2)" x R^™. The Lagrangian is 

L{9, g, 1, i) = a^(7) + + Tb. 

The Lagrangian here is partly degenerate, i.e., it is affine in point vortex velocities. Thus, 
for numerical purposes, it fits into th e framework of variational integrators for degenerate 



Lagrangian systems, see iRowlev and Marsden. (,2002. ). 



5 Conclusions and Outlook 

We have introduced the Hamiltonian description for several rigid bodies interacting with 
point vortices, assuming zero circulation around the individual bodies, but arbitrary point 
vortex strengths. This system has been known only for the case of a single rigid body. 
We have used the general framework of cotangent bundle reduction only to determine the 
reduced phase space of the system, as well as the general structure of the symplectic form 
on the reduced phase space. From there we have determined the symplectic form directly, 
without resorting to the abstract framework of mechanical connections. We have derived 
the equations of motion for the general case, and also for the Lie-Poisson case of a single 
rigid body with zero total strength of the point vortices. From the Hamiltonian formula- 
tion we have also given a Lagrangia n description of the dynam ics. This allows to use the 



framework of variational integration (jMarsden and Westl . 120011 ) to construct structure- and 
momentum-preserving integrators. 

We expect that our formulation generalises to the 3-dimensional case, describing the 
dynamics of several rigid bodies interacti ng with vortex filaments. This system is also only 
known for the case of a single rigid body ( Shashikanth et al. . 20081) . 



Acknowledgment: It is a pleasure to thank Ulrich Pinkall for his support, for inter- 
esting discussions and invaluable suggestions. This work is supported by the DFG Research 
Center Matheon and the SFB/TR 109 "Discretization in Geometry and Dynamics". 

References 

Arnold, V. (1966). Sur la gomtrie diffrentielle des groupes de lie de dimension infinie et 

ses applications I'hydrodynamique des fluides parfaits. Annates de I'institut Fourier, 
16(1):319-361. 



10 



Borisov, A. V., Mamacv, I. S., and Ramodanov, S. M. (2007). Dynamic interaction of point 
vortices and a two-dimensional cylinder. J. Math. Phys., 48(6):065403. 

Kirchhoff, G. R. (1870). Uber die Bewegung eines Rotationskorpers in einer Fliissigkeit. 
Journal fiir die Reine und Angewandte Mathematik, 71:237-262. 

Lin, C. C. (1941). On the motion of vortices in two dimensions - I and II. Proc. Natl. Acad. 
Sci. U.S.A., 27:570-575. 

Marsden, J., Misiolek, G., and Ortega, J. (2007). Hamiltonian Reduction by Stages. Lecture 
Notes in Mathematics. Springer, Berlin. 

Marsden, J. E. and West, M. (2001). Discrete mechanics and variational integrators. Acta 
Numer., 10:357-514. 

Nair, S. and Kanso, E. (2007). Hydrodynamically coupled rigid bodies. J. Fluid Mech., 
592:393 411. 

Rowley, C. W. and Marsden, J. E. (2002). Variational integrators for degenerate La- 
grangians, with application to point vortices. Proc. 41st IEEE Conference on Decision 
and Control, 2:1521-1527. 

Saffman, P. G. (1992). Vortex Dynamics. Cambridge Univ. Press. 

Shashikanth, B. N. (2005). Poisson brackets for the dynamically interacting system of a 2D 
rigid cylinder and N point vortices: The case of arbitrary smooth cylinder shapes. Regul. 

Chaotic Dyn., 10:1 14. 

Shashikanth, B. N., Marsden, J. E., Burdick, J. W., and Kelly, S. D. (2002). The Hamiltonian 
structure of a two-dimensional rigid circular cylinder interacting dynamically with N point 
vortices. Phys. Fluids, 14(3):1214. 

Shashikanth, B. N., Sheshmani, A., David, S., and Jcrrold, K. (2008). Hamiltonian structure 
for a neutrally buoyant rigid body interacting with N vortex rings of arbitrary shape: the 
case of arbitrary smooth body shape. Theor. Comput. Fluid Dyn., pages 37-64. 

Vankcrschaver, J., Kanso, E., and Marsden, J. E. (2009). The geometry and dynamics of 
interacting rigid bodies and point vortices. J Geom Mech., l(2):223-266. 



11 



